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The geometric phase can act as a signature for critical regions of interacting spin chains in the limit where the 
corresponding circuit in parameter space is shrunk to a point and the number of spins is extended to infinity; for 
finite circuit radii or finite spin chain lengths, the geometric phase is always trivial (a multiple of 2n). In this 
work, by contrast, two related signatures of criticality are proposed which obey finite-size scaling and which 
circumvent the need for assuming any unphysical limits. They are based on the notion of the Bargmann invariant 
whose phase may be regarded as a discretized version of Berry's phase. As circuits are considered which are 
composed of a discrete, finite set of vertices in parameter space, they are able to pass directly through a critical 
point, rather than having to circumnavigate it. The proposed mechanism is shown to provide a diagnostic tool 
for criticality in the case of a given non-solvable one-dimensional spin chain with nearest-neighbour interactions 
in the presence of an external magnetic field. 

PACS numbers: 03.67.Mn, 05.70.-a 



I. INTRODUCTION 
A. Quantum Phase Transitions 

Unlike ordinary phase transitions, quantum phase transi- 
tions (QPTs) 0]] take place at a temperature of absolute zero 
and are consequently driven entirely by quantum fluctuations, 
rather than thermal ones. Moreover, the long-range, alge- 
braically decaying correlations which are characteristic of 
critical many-body ground states are due entirely to entan- 
glement. As such, one would expect entanglement measures 
(see Ref. JUl for an overview) to provide further insight into 
the fundamental physical underpinnings of QPTs, possibly 
complementing the conventional condensed-matter approach, 
which relies mostly on two-point correlation functions. 

Indeed, Osterloh and co-workers JH] demonstrated singular 
and scaling behaviour of a bipartite entanglement measure in 
the vicinity of the critical point in the one-dimensional XY 
spin model. On the other hand, this behaviour could be under- 
stood entirely from that of the relevant two-point correlation 
functions as these are sufficient to determine the two-particle 
reduced density matrices. Furthermore, the entanglement be- 
tween just two sites is not particularly well suited to capturing 
the large-scale behaviour of correlations, which becomes all 
the more relevant in the critical regime. This was recognized 
in Ref. J4[], where the scaling behaviour of the entanglement 
between a contiguous block of L sites with the rest of the lat- 
tice was first considered in lattice field theories. It was found 
that the entanglement entropy of noncritical chains with short- 
ranged interactions reaches a saturation level, while the entan- 
glement entropy diverges logarithmically with L in the field 
limit, i.e. in the case of critical chains. This approach was 
subsequently pursued numerically in the particular context 
of the one-dimensional XY model and, in more generality, in 
an analytical treatment based on random matrix theory @|, 
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applicable to any spin chain Hamiltonian that can be cast into 
a quadratic form of fermionic operators. 

At the same time, progress has been hampered by the re- 
alization that bi-partite entanglement, whether of blocks of 
spins or otherwise, at best only gives us an incomplete, lo- 
cal picture of the entanglement exhibited by generic many- 
body ground states. It is something of a truism, of course, to 
note that a truly global, multi-partite approach is indispens- 
able if one wishes to capture the entanglement that pervades 
critical many-body ground states simultaneously at all length 
scales (for recent work on multi -partite entanglement in quan- 
tum spin chains see for example Refs. 10, [Hi)- 

Unfortunately, the theory of multi-partite entanglement is 
still very much in its infancy. It is true that many of the en- 
tanglement measures used for bi-partite states carry straight- 
forward generalizations to the multi-partite setting. This is 
particularly true for distance-based entanglement measures, 
such as the relative entropy of en tang lement ^ or the closely 
related geometric entanglement fllQfl : the entanglement of a 
state is then simply quantified in terms of the minimum dis- 
tance of that state from the set of all multi-partite separable 
states, rather than the set of all bi-partite separable states. 
However, it is one thing to generalize the axiomatic defini- 
tions of entanglement measures to include the multi-partite 
case, but quite another to still be able to compute these in 
practice. Moreover, the theory of multi-partite entanglement 
is still plagued by a host of different candidates for suitable 
entanglement measures, even for the case of pure states (see 
Ref. yfl for more details). 

In view of these difficulties, attention has shifted to include 
other, potentially related, means of characterizing QPTs ifTHl . 
One such approach centres around the notion of geometric 
phase d El BEES, 

and it is this approach which we 
wish to pursue in the remainder of this work. 
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B. Geometric Phase 



In his seminal paper 11211 . Berry investigated the phase 
picked up by an eigenstate of a parameter-dependent Hamilto- 
nian when transported adiabatically around a closed trajectory 
in parameter space. It turns out that in addition to the well- 
known dynamical phase there is also a geometric component 
to the phase, which is observable, at least in principle. While 
the dynamical phase provides a measure of the duration of the 
Hamiltonian's evolution and is independent of the geometry 
of the trajectory followed, conversely, the geometric phase is 
independent of the rate at which the state is transported around 
the loop (as long as this is slow enough for the adiabatic theo- 
rem to rule out transitions to neighbouring, orthogonal states) 
and depends solely upon the geometry of the trajectory. The 
geometric phase had hitherto been widely overlooked as just 
another unphysical phase factor, and certainly had not been 
granted the level of recognition it enjoys nowadays. 

Before going on to explain how Berry's phase may be used 
to probe for QPTs, it is useful here to first restate a simple 
example given in Berry's original paper, which serves to illus- 
trate the concept of geometric phase. To that end, consider a 
single spin- 1/2 particle which is coupled to an external mag- 
netic field, and suppose the spin is initially aligned with the 
magnetic field, i.e. the particle is in an eigenstate of the sys- 
tem's Hamiltonian. Then the spin will remain aligned with 
the magnetic field vector (the particle remains in the instan- 
taneous eigenstate) when the magnetic field vector is made to 
rotate adiabatically. Upon completion of some closed trajec- 
tory, the final state of the particle differs from its initial one 
by an overall phase factor, which is part dynamical, part ge- 
ometric in origin. Now, it turns out that the geometric phase 
is in fact equal in size to precisely half the solid angle sub- 
tended at the origin by the magnetic field vector's trajectory 
in parameter space. This result not only serves to highlight 
in a particularly acute way the geometric character of Berry's 
phase, but it also gives an indication of the origin of Berry's 
phase: the spin's state is degenerate when the magnetic field is 
switched off, and the geometric phase can thus be interpreted 
as providing a measure of the view of the circuit as seen from 
that point of degeneracy at the origin of the parameter space. 

In fact, in any general setting, the emergence of geometric 
phases can always be traced back to the presence of isolated 
singularities in parameter space. Formally, Berry's phase can 
be related to the curvature of the Hilbert space bundle over the 
base space of parameters fl7il . and it is this interpretation of 
Berry's phase which renders its potential usefulness as a diag- 
nostic tool for criticality plausible: points of degeneracy are 
associated with a greater curvature of the associated Hilbert 
space bundle, and therefore it might be reasonable to expect 
Berry's phase to act as a signature for quantum critical points 
in interacting spin systems. 

Of course, the geometric phase of N non-interacting spin- 
1/2 particles with respect to a particular circuit is simply N 
times the geometric phase of a single particle, which in turn 
can be related to the circuit's solid angle, as explained above. 
However, as we switch on the spin-spin interactions, this sim- 
ple interpretation of the geometric phase in terms of solid an- 



gles starts to break down. It is of great interest to know how 
the geometric phase of a set of interacting spins with respect 
to a given circuit changes as a function of the coupling pa- 
rameter, and whether the geometric phase is able to signal 
the presence of critical points in the system. An important 
point to note here is that the circuit in parameter space need 
only pass near the critical point for the Berry phase to regis- 
ter it. In other words, the system need not actually undergo 
the quantum phase transition for the Berry phase to pinpoint 
its presence and location in parameter space, a consideration 
which assumes particular importance in the light of the diffi- 
culties that may be associated with physically implementing 
actual QPTs. 



C. Previous Work 

A relationship between the geometric phase and criticality 
in spin chains was examined by Carollo and Pachos lfl3i[l4ll . 
Specifically, the authors analyzed the Berry phase of the 
ground state of a one-dimensional spin- 1/2 XY model 
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where M = (N — l)/2 for an odd number of spins N. The 
Berry phase, as computed with respect to a rotation of the 
complete Hamiltonian around the z-axis, was shown to exhibit 
the following property in the thermodynamical limit N — > oo: 
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As the XY model is critical for |A| < 1 (and non-critical 
otherwise), the quantity given above may be interpreted as 
a signature of criticality for this model. Unfortunately, 
the Berry phase per spin, <p/M, is not a physical quan- 
tity that is accessible via experiments. We can only mea- 
sure the total Berry phase (p, which vanishes for all val- 
ues of A: i.e. (1/M) lim 7 ^o <£. = (mod 2ir), VA. This 
is in agreement with Hamma 1151 . who also showed, via 
a different line of argument, that for a finite number of 
spins, the Berry phase acquired by the ground state is al- 
ways trivial (0 or 2tt). As an alternative, it was suggested 
that the same quantity would be non-trivial in the thermo- 
dynamic limit: lim 7 ^o \vn\M^oo{f/M) 7^ 0. However, 
for any physical, finite system, no matter how large, one is 
not able to detect criticality in this manner. In other words, 
any finite-size scaling behaviour is completely lacking, i.e. 
limM-»oo YmLy^fy / M) = (mod 2tt); it is only in the ther- 
modynamic limit itself that this method represents a signature 
of criticality, i.e. the limiting case is reached discontinuously. 
As such, the result's worth is perhaps more of an abstract, 
mathematical nature, than of any real, physically measurable 
consequence. 

In the present work, in contrast, we employ a technique 
based on the geometric phase that is able to detect critical 
points without the need to contract the circuit's radius to zero 
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or to extend the number of particles to infinity. We do so by 
considering discretized circuits that pass directly through the 
critical point, rather than circumnavigating it. This procedure 
and the obtained results are outlined in the following. 



tices of the circuit. Of particular interest, of course, is what, 
if anything, may be said about the Bargmann invariant as the 
coupling parameter crosses its critical point J c . The report 
concludes with a summary and brief discussion of our work. 



H. PROCEDURE AND OUTLINE 

We start by introducing the concept of a Bargmann invari- 
ant jlHl and its associated phase, which may be regarded as 
a generalized Berry phase. It will be seen that Bargmann 
invariants are ideally suited for numerically analyzing the 
gauge-invariant phase induced by evolving states along a dis- 
cretized circuit in parameter space. We then introduce the 
model of interacting spin- 1/2 particles with which we will 
be concerned for the rest of this paper and discuss some of 
its more salient features. The Bargmann invariant and associ- 
ated phase is computed numerically for a chosen discretized 
circuit and plotted as a function of the spin coupling strength 
J. Of note will be the 'speed' (to be defined) with which the 
Bargmann invariant changes as a function of J and the be- 
haviour of the magnitude of the Bargmann invariant; we will 
find that both of these quantities can act as adequate signa- 
tures of criticality. At this stage, it may be helpful to depict the 
procedure schematically: Let us fix the coupling parameter to 




FIG. 1 : A schematic representation of the procedure used to compute 
the Bargmann invariant with respect to some chosen circuit {g_ B } in 
parameter space, as a function of the spin coupling strength J. 

some value Jo and choose a particular circuit composed of JV 
vertices a s , each representing the set of parameters needed 
(together with Jo) in order to completely describe the Hamil- 
tonian of the system at that point. At each individual point 
on the circuit we numerically compute the ground state of the 
system. The Bargmann invariant associated with the circuit at 
Jo is then obtained by cyclically multiplying successive over- 
laps of these ground state wavefunctions. This will be ex- 
plained in more detail in section|III] which is devoted entirely 
to the subject of Bargmann invariants; for now, we merely 
note that, in order to be sure to acquire a gauge-invariant quan- 
tity in this way, the wavefunctions at the start and end points of 
the circle need to be taken as identical. By repeating this pro- 
cedure for a range of different coupling strengths (generating 
a complex number, the Bargmann invariant, for each value), 
we start to build up a picture of the general trend. We present 
our results for varying lengths of the spin chain, contrast the 
case of spin- 1/2 particles with that of spin-1 particles, and dis- 
cuss the effect of increasing the number N of constituent ver- 



m. BARGMANN INVARIANT AND PHASE 

In quite general terms one may define the phase difference 
X between two non-orthogonal state vectors, and \^>2), 
by the relation 

cxp(ix) = (*l|* 2 )/|(*l|*2)|, 

so that x = aigi^ i\^ 2) ■ Naturally, individual state vectors 
are only defined up to arbitrary phase factors; x 1S thus gauge- 
dependent and its value may not be assigned any direct phys- 
ical meaning. On the other hand, any cyclic combination of 
state vectors, such as arg((^I>i | ^2 } { , 3 / 2 1 ^3) (^3 i'S'i)) is mani- 
festly gauge-invariant and may therefore potentially represent 
a physically relevant, measurable quantity. This observation 
features strongly in the early works of Pancharatnam lfl9ll and 
it also appears in the form of a remark in the celebrated proof 
of Wigner's theorem by Bargmann ifTill . 

Introducing some notation, we define the Bargmann invari- 
ant with respect to a general TV-vertex circuit in parameter 
space by 

JV-1 AT-l 
C = n (H^s+i) = tr J] \A)&s\, i>M = 1*. (2) 

s=0 s=0 

The associated Bargmann phase is denoted by ip = arg(C). 
Modulo 2ir, the global Bargmann phase is just the sum of the 
individual phases, i.e. if = X^o* ar g(V's|V's+i)- Note that 
we require M 3 for non-trivial Bargmann phases. 

Suppose now, for argument's sake, that we are dealing 
with state vectors \i[)(a s )) that represent the instantaneous 
ground states of a family of parameter-dependent Hamiltoni- 
ans H({a}). Then, in the continuum limit where the sum 
above is replaced by a contour integral over an infinitesimal 
phase along a (smooth) closed circuit in parameter space, the 
Bargmann invariant reduces to the usual Berry phase. In this 
sense Bargmann invariants may be regarded as generalized 
Berry phases. The key advantage afforded by the generalized 
formulation lies in its computational ease: Bargmann invari- 
ants lend themselves directly to numerical computation, and 
thus, crucially, to scenarios where parts of the circuit under 
consideration represent non-integrable Hamiltonians, as is the 
case in the present work. An obvious potential drawback is 
that the discretization procedure may lack the simple interpre- 
tational underpinnings of the Berry phase in terms of physical 
adiabatic processes. 

IV. THE SPIN MODEL 

Consider a spin chain with nearest-neighbour XX inter- 
actions in the presence of an external magnetic field, which 
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is aligned along an arbitrary direction n . The system is de- 
scribed by the following Hamiltonian: 

N N 

H J t = jJ2°M+i+Bj2<rk, N>2, (3) 

k=l k=l 

where, n = (n x ,n y ,n z ) denotes a vector of unit length, so 
that a " = J2 a n a° a with J2 a n l = 1, a € {x, y, z}. Note 
that we impose periodic boundary conditions, so that oj^+i = 
of. 

As already outlined in the previous two sections, our aim is 
now to analyze the Bargmann invariant induced by the ground 
states of a family of Hamiltonians (0 that is characterized by 
a set of N unit vectors {it}. Typically, we will imagine the 
magnetic field vectors of that family of Hamiltonians to be 
arranged on the vertices of a regular polygon residing on the 
surface of the Bloch sphere (i.e. B = 1); an example of such a 
circuit is depicted in Fig. [2] For a given coupling strength, Jo, 
the Bargmann invariant can now be obtained by numerically 
computing the ground states at each of the circuit's vertices 
and applying Eq. (0. Having chosen our circuit in parameter 
space, we would then like to investigate the manner in which 
the Bargmann invariant depends on the choice of spin cou- 
pling strength, J. Of particular interest are circuits that slice 
through any critical points of the system. Does the Bargmann 
invariant undergo a marked shift when the parameters pass 
through their critical values? Of further interest is the depen- 
dence on the number N of vertices of the regular polygon 
circuit, as well as the dependence on the number N of spins 
in the chain. 

Leaving the choice of circuit in parameter space open for 
now, it is evident from the form of Hamiltonian (01 that there 
are two 'special' directions in which the magnetic field vector 
may point, namely the x-direction and any direction in the y — 
z plane. The former is essentially a classical model while the 
latter corresponds to the quantum Ising model. In both cases, 
the Hamiltonian is analytically soluble. Of course, it is well 
known that the Ising model exhibits criticality when B = J, 
provided we find ourselves in the thermodynamic limit of an 
infinite spin chain. Aside from the two special cases outlined, 
however, the Hamiltonian will in general be non-integrable at 
any point along the circuit; this forces us to resort to numerical 
simulations of a limited number of spins, which moves the 
Ising model's critical point outside our region of accessibility. 
On the other hand, when the magnetic field vector points in 
the x-direction, the Hamiltonian does possess a critical point, 
even for a finite number of spins, as outlined in the following. 

V. THE CRITICAL POINT 



following we will demonstrate that H x possesses a critical 
point JlOD at J = J c = \B\/2. 

It is not difficult to find the ground state of this model. 
When J ^ J c , it is simply 

(0) (\-f N , forB^0&J< J c , 
1 ' \\+)® N , forB<0&J< J c , 

with the corresponding (non-degenerate) ground state energy 
given by = N(J - 2J C ). 

In the regime J ^ J c , the situation is only slightly more 
involved as the parity of N now becomes important. For an 
even number of spins, the (doubly degenerate) ground state 
energy is given by E^ = —NJ, the corresponding eigen- 
states being | H — ^® N / 2 an( j | — ^®n/2 ^ p Qr an o ^ nurm -, er 
of spins, the ground state energy is A^-fold degenerate and 
reads E^ Q > = —NJ+2(J — J c ). The ground state is given by 




- + - + +-}, for B ^ & J ^ J c , 

+ - + - + ••• + -+), for B < & J ^ J c , 



and all N translations thereof. In order to convince oneself 
that this is indeed the correct ground state, one can easily show 
that no individual spin flip is capable of lowering the energy 
any further. 

From the analysis above it follows that the ground state en- 
ergy £(°> can always be written as a function of J — J c |, 
demonstrating non-analyticity at the critical point J c . More- 
over, the energy gap vanishes in the region J ^ J c . 

VI. THE CIRCUIT 

In light of the previous discussion, the parameter circuit 
which suggests itself is the following: a regular polygon re- 
siding on the Bloch sphere of which one edge is bisected by 
the x-axis, as illustrated in Fig. [2] As we move along the ver- 




A special case of the class of Hamiltonians ® occurs when 
the magnetic field vector points in the x-direction: 



N 



N 



H x = jJ2^k + i + Bj2^k, N^2. 



(4) 



k=l 



This Hamiltonian is purely classical: its two constituent terms 
commute and are diagonal in the {|±}} eigenbasis. In the 



FIG. 2: A schematic representation of the circuit's topology. The 
critical x-axis slices mid-way through an edge that connects two ad- 
jacent vertices of the polygon. 
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tices of the circuit, the critical x-axis is crossed in the direction 
of positive z-axis. 



VII. RESULTS FOR THE SPIN-1 /2 CHAIN 

We now turn to the results of our simulation of the 
Bargmann invariant for the spin-1/2 chain, which are sum- 
marized in Fig. [3] The circuit used in the simulation is the 
regular polygon described in the previous section, but con- 
sists of N = 100 vertices. The Bargmann invariant and its 
phase were computed for different spin chain lengths, ranging 
in odd numbers from 3 to 11, and shown in Fig.|3]as increas- 
ing from the top. Note that the results for even numbers of 
spins have been omitted as they showed no qualitative differ- 
ence. This may be due to finite-size effects of the kind which 
have been observed in similar studies elsewhere, see for exam- 
ple Ref . ll2lll . The left-hand column of the figure depicts the 
Bargmann invariants C on the unit disk in the complex plane. 
Each data point of the graph stems from a different value of 
the coupling parameter J. The right-hand column shows the 
corresponding Bargmann phase ip/ir, as a function of the cou- 
pling parameter (J — J c ) / r, where r refers to the radius of the 
circle within which the TV-polygon circuit is inscribed. Note 
that in this case the radius was chosen as 10~ 5 , but the graphs, 
as they are plotted, are qualitatively independent of the size of 
this radius. 

We would like to point out several features that are appar- 
ent from Fig. [3] First and foremost, the Bargmann invariant C 
is observed to always pass through the origin of the complex 
plane just when the spin coupling strength attains its critical 
value of J c = 1/2. In other words, the magnitude \C\ van- 
ishes at the critical point. This is one of the features which 
we would like to propose as a signature of criticality, and is 
discussed further in section lVIII Al 

Another trend that is immediately apparent from Fig. [3] is 
that with increasing spin chain length the Bargamann invari- 
ants 'wind ever more loops' around the origin. This statement 
can be made a little more precise by referring to the corre- 
sponding Bargmann phases. These graphs are characterized 
by discontinuities at the critical point: 'jumps' by it alternat- 
ing with 'jumps' by 2ir. However, on joining up those separate 
pieces in the graphs end-to-end, one notices that the 'image' 
covered by the resulting graphs grows with N as ir(N — l)/2. 
This result may be explained by way of referring to the J = 
limit of non-interacting spins for which the Berry phase grows 
linearly with the number of spins. Also note that for spin-1 
chains, the 'joined-up' Bargmann phase covers twice as much 
ground for any given number of particles (see Fig.|7J. 

In addition, we should note that we obtain Berry's solid 
angle result for the case of non-interacting spins, J — 0. Ap- 
proximating the solid angle of our regular polygon circuit as 
that of the enveloping cone, the Bargmann phase ip asymptot- 
ically approaches Berry's solid angle result as the number of 
vertices of the polygon is increased. 

Of course, in the present work only spin chains of up to 11 
particles were simulated, and the exponential growth in the 
size of the Hilbert space with increasing particle number pre- 
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FIG. 3: Bargmann invariants and their phases. The rows corre- 
spond to the number of spin-i particles N, increasing from the top as 
3, 5, 7, 9, 11. The left column depicts the Bargmann invariants C on 
the unit disk in the complex plane, each point of the graph stemming 
from a different value of the coupling parameter J. The right-hand 
column shows the corresponding Bargmann phase tp/n as a function 
of the coupling parameter (J — Jc)/r. See the text for a detailed 
description of observations. 



eludes one from simulating chains that are much longer than 
that I22I1 . However, from the preceding discussion the trend 
for Bargmann invariants and phases of longer spin chains 
should by now be fairly obvious. In particular, it has be- 
come clear that, perhaps surprisingly, precious little may be 
deduced about criticality from the Bargmann phase. Instead, 
a much clearer indication of the existence of a critical point is 
the vanishing magnitude of the Bargmann invariant itself. 
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VIII. SIGNATURES OF CRITICALITY 
A. Magnitude of the Bargmann Invariant 

It was already apparent from Fig. [3] that the magnitude of 
the Bargmann invariant C is able to assume the role of a diag- 
nostic tool for critical points in the finite-sized interacting spin 
chain considered. For clarity, the magnitude of the graphs of 
C in Fig. [5] are plotted separately in Fig.|U An analogous plot 




-2 -1.5 -1 -0.5 0.5 1 1.5 



(J - J c )/r 

FIG. 4: Plots of \C\ for varying lengths of the spin chain, ranging 
from 7V = 3to7V = ll (odd numbers only). The longer the spin 
chain, the more sharply peaked is the graph. The number of vertices 
on the circuit is kept constant, at N — 100. 

is shown in Fig. [5] but this time for a fixed number of spins N 
and a varying number of circuit vertices W. Of course, there 




0.2 
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FIG. 5: Plots of \C\ for varying numbers of vertices on the circuit, 
ranging from N = 100 to N = 300 (in steps of 50). The more 
vertices there are on the circuit, the more sharply peaked is the graph. 
The spin chain length is kept constant, at just N = 3. Note that 
qualitatively nothing changes as we increase the chain length (tested 
up to TV = 13). 

is a perfectly plausible explanation for these results: as the 
circuit passes over the critical point the ground state changes 
abruptly, so that the overlap of the ground states on either side 
of the critical point will be rather small in magnitude. This 
in turn forces the magnitude of the product of overlaps along 



the circuit in Eq. [2] to be small. In fact, similar results were 
observed in a work by Zanardi and Paunkovic ll ill , who con- 
cluded that the ground state overlap function is itself a good 
characterization of QPTs. 

B. 'Speed' of the Bargmann Invariant 

Another characteristic that is sensitive to the phase transi- 
tion is the 'speed' with which the Bargmann invariant changes 
as a function of the coupling parameter, which we define by 
v s := \C a +i — C s \/(J s +i — Js)- As is evident from Fig. [6] the 
speed picks up markedly in the vicinity of the critical point. 
Again, this phenomenon is a direct result of the rapid change 
of the ground state near a critical point. 
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FIG. 6: Plot of the 'speed', v s := \C a +i — C s \/(J s +i — Js), as 
a function of the spin coupling strength. The spin chain length is 
TV = 3 and the circuit is composed of N = 100 vertices. 



IX. RESULTS FOR THE SPIN-1 CHAIN 

The results presented in Fig.|7]derive from simulations that 
are entirely analogous to those that gave rise to Fig. [3] but 
they refer to the spin-1 chain model, rather than the spin- 1/2 
model, Eq. [3] The Bargmann invariants are shown for spin- 
1 chains of lengths N — 5 (left-hand column) and N = 7 
(right-hand column), and allow for direct comparison and con- 
trast with the spin- 1/2 case of Fig. [3] The circuit is composed 
of M = 100 vertices. Again, the uppermost row shows the 
trajectory traced out in the complex plane by the Bargmann 
invariant as the spin coupling strength is varied from J c — 2r 
to J c + 2r. The second row from the top shows the corre- 
sponding Bargmann phase; notice that points in the immedi- 
ate vicinity of the critical point have been disregarded as the 
phase of a complex number becomes increasingly prone to 
error as the origin of the complex plane is approached. The 
next row pictures the overall 'extent' of the phase, with the 
individual pieces of the graph having been joined up end-to- 
end, as described in section [VTT1 Finally, the magnitude of the 
Bargmann invariant is shown, and is clearly seen to vanish in 
the vicinity of the critical point. 
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FIG. 7: These plots depict the Bargmann invariants for spin-1 chains 
of lengths N = 5 (left-hand column) and N = 7 (right-hand col- 
umn), and allow for direct comparison and contrast with the spin- 1/2 
case of Fig. [5] For a detailed description of the results the reader is 
kindly referred to the text in section Hxl 



Qualitatively, the results for spin-1 chains are very simi- 
lar to those of spin- 1/2 chains. In particular, the magnitude 
of the Bargmann invariant still represents a signature for the 
critical point. Of note is also the fact that the 'joined-up' 
Bargmann phase extends twice as high for spin-1 chains as 



it does for spin-1/2 chains, for the same number of particles. 
Again, this can be understood from Berry's solid angle re- 
sult for non-interacting particles, where the geometric phase 
is proportional to the spin dimension. 



X. SUMMARY AND DISCUSSION 

In summary, we have shown how the magnitude and 'speed' 
of the Bargmann invariant can act as a signature of the critical 
point for the finite one-dimensional spin chain model consid- 
ered. This stands in stark contrast to the Bargmann phase, 
which is oblivious to the phase transition as long as the spin 
chain length remains finite. 

It is conjectured that the signatures presented here in the 
context of a specific spin chain model will uphold their valid- 
ity also in a wider, more general setting. It would be worth- 
while to conduct experimental investigations with a view to 
testing this conjecture; similarly, any theoretical proof of the 
proposed conjecture would be of considerable interest. 

Finally, it is worth stressing that the observed signatures 
of criticality are very likely a direct consequence of the finite 
'speed' with which the circuit is traversed in parameter space, 
i.e. a result of the discrete nature of the circuit. No matter how 
finely spaced are the neighbouring vertices on the circuit, the 
adiabatic approximation will always break down in the im- 
mediate vicinity of the critical point, causing Landau-Zener 
tunnelling effects ll23ll to take centre stage. In the limit of a 
smooth, continuous circuit, the observed signatures of criti- 
cality would likely disappear altogether. 
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